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The study of the Sivers effect, describing correlations between the transverse polarization of 
the nucleon and its constituent (unpolarized) parton’s transverse momentum, has been the topic 
of a great deal of experimental, phenomenological and theoretical effort in recent years. Semi- 
inclusive deep inelastic scattering measurements of the corresponding single spin asymmetries (SSA) 
at the upcoming CL AS 12 experiment at JLab and the proposed Electron-Ion Collider will help 
to pinpoint the flavor structure and the momentum dependence of the Sivers parton distribution 
function describing this effect. Here we describe a modified version of the PYTHIA6 Monte Carlo 
event generator that includes the Sivers effect. Then we use it to estimate the size of these SSAs, in 
the kinematics of these experiments, for both one and two hadron final states of pions and kaons. 

For this purpose we utilize the existing Sivers parton distribution function (PDF) parametrization 
extracted from HERMES and COMPASS experiments. Using this modified version of PYTHIA6, we 
also show that the the leading order approximation commonly used in such extractions may provide 
significantly underestimated values of Sivers PDFs, as in our Monte Carlo simulations the omitted 
parton showers and non-DIS processes play an important role in these SSAs, for example in the 
COMPASS kinematics. 

PACS numbers: 13.88.+e, 13.60.-r, 13.60.Hb, 13.60.Le 

I. INTRODUCTION 

The study of the structure of the nucleon has long been 
a focus of both theory and experiment within the medium 
and the high energy hadronic physics community, mostly 
concentrating on the distributions of the longitudinal mo¬ 
mentum fraction of partons within nucleon described by 
the transverse momentum integrated (’’collinear”) par- 
ton distribution functions (PDF). Lately, exploring the 
transverse structure of the nucleons in momentum space 
via transverse momentum dependent (TMD) PDFs and 
the generalized parton distribution functions in impact 
parameter space has attracted great interest. The spin 
of the nucleon and its constituent partons plays an impor¬ 
tant role when describing the transverse structure of the 
nucleon. Here we will focus on the so-called Sivers TMD 
PDF that describes the modification of the fully unpo¬ 
larized TMD PDF when considering unpolarized partons 
inside of a transversely polarized nucleon [I] . The knowl¬ 
edge of the Sivers PDF is crucial in our understanding of 
the spin structure of the nucleon. It can be accessed in a 
variety of experiments in the so-called deep inelastic scat¬ 
tering (DIS) regime, where the probes intact directly with 
a parton inside of the target, such as the Drell-Yan pro¬ 
cess and semi-inclusive hadron production in polarized 
hadron-hadron collisions, or semi-inclusive DIS (SIDIS) 
of a lepton off a transversely polarized nucleon. 


In a SIDIS process with a single detected final state 
hadron, the Sivers PDF can be accessed through a term 
in the single spin asymmetry (SSA) modulated with the 
sine of the so-called Sivers angle [2], defined as the az¬ 
imuthal angle between the transverse momentum of the 
detected hadron and the transverse spin of the nucleon. 
The amplitude of this modulation is given as a convo¬ 
lution of the Sivers PDF with the TMD unpolarized 
fragmentation function (FF) divided by a convolution 
of the unpolarized TMD PDF and TMD FF. Experi¬ 
mental measurements of these asymmetries have been 
performed by several collaborations: HERMES mm, 
COMPASS [SIS] and JLab HALL A [7], extracting siz¬ 
able effects for both pion and kaon final states. The 
Sivers PDF has then been extracted in Ref. [Sj using 
phenomenological parametrizations of the involved TMD 
PDFs and TMD FFs to fit the experimentally measured 
SSAs, and more recently with the inclusion of the TMD 
QCD evolution 0 in Refs. mm- Unfortunately these 
parametrizations suffer from relatively large uncertain¬ 
ties arising from both theory and experiment. These 
uncertainties come from the experimental measurements 
themselves, the lack of knowledge of unpolarized TMD 
FFs (especially those for the unfavored channels [12], as 
well as those involving the strange quark [13| ). and as¬ 
sumptions about the shape and the flavor dependences 
of both unpolarized and Sivers PDFs and unpolarized 
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FFs. This is relevant, as in calculating SSAs one inte¬ 
grates over the transverse momentum and then sums over 
the type of struck parton. Moreover, in such extractions 
there are unquantified uncertainties associated with us¬ 
ing the leading order expressions for the SSAs that ignore 
any parton showering (note that TMD evolution partially 
accounts for this effect) or non-DIS process contributions 
(such as vector-meson-dominance) to the cross section of 
the inclusive hadron production in lepton scattering off 
a polarized nucleon target. 

Recently, it has been suggested that the Sivers PDF 
can be also measured in SIDIS production of two hadrons 
in the fragmentation of the same struck parton mm- 
Here the SSA has two terms involving the Sivers PDF, 
modulated with respect to the sine of the Sivers angles 
of the total and the relative transverse momenta of the 
pair. In each of these terms, the Sivers PDF is con¬ 
voluted with a corresponding unpolarized dihadron frag¬ 
mentation function (DiFF), that is not yet known. LEPTO 
Monte Carlo (MC) generator was modified to study these 
SSAs by using the Sivers PDF parametrizations of pun] 
and the Lund string fragmentation model employed in 
LEPTO to describe the semi-inclusive production of two 
hadrons. Using this MC generator, it was shown that 
the corresponding SSAs are of a similar size to those for 
the single hadron production in the kinematics of the 
COMPASS experiment. In addition, the two hadron fi¬ 
nal state provides a larger basis than the single hadron 
case for disentangling the flavor dependence of the Sivers 
PDF, since for N types of detected hadrons we can mea¬ 
sure N single hadron SSAs compared with N ■ (N + l)/2 
dihadron SSAs (including those with the same type of 
hadrons in the pair). Typically, due to the design con¬ 
straints of a particular detector, N is not large enough 
to determine the flavor dependences. On the other hand, 
as shown in Sec. [Ill] and Sec. |IV| a particular choice of 
hadrons in a pair may allow one to gain more sensitivity 
to the contributions of a quark with given flavor. 

Thus the measured dihadron SSAs will help to put 
additional constraints on the flavor dependence of the 
Sivers PDFs, and ideally extract them from their convo¬ 
luted sums in SSAs without any assumptions on their fla¬ 
vor dependences. Moreover, SSA measurements in two- 
hadron SIDIS will be actively pursued in the upcoming 
CLAS12 and proposed SoLID experiments for accessing 
the transversity PDF through the interference fragmen¬ 
tation channel Ellis, making Sivers two-hadron studies 
a relatively easy addition to the analysis program. To¬ 
gether with the one hadron SIDIS and other measure¬ 
ments this will help to map out the full momentum and 
flavor dependence of Sivers PDFs. This is, of course, crit¬ 
ically hinged on our detailed knowledge of unpolarized 
TMD PDFs, TMD FFs and DiFFs, that also need to 
be mapped out in future experiments. Finally, dihadron 
measurements will help to investigate the correlations in 
the hadronization process that have been discussed re¬ 
cently [B5J. 

The exploration of the transverse structure of the nu¬ 


cleon, and the Sivers function in particular, is one of the 
main focuses of the CLAS12 experiment at the soon to be 
operational JLab at 12 GeV and the proposed Electron 
Ion Collider (EIC) [2D]. Thus it is important to make 
projections for the sizes of the expected SSAs using the 
current information on the Sivers PDF extracted from 
experiment, as well as to analyze what are the important 
physics questions that can be posed in different centre-of- 
mass energy regimes. Here, we use a modified version of 
the PYTHIA6 (mPYTHIA) MC event generato r th at incor¬ 
porates the Sivers effect (as described in Sec. Ill) to make 
such projections for the CLAS12 experiment at JLab and 
for two sets of beam energies for EIC. 

This paper is organized in the following way. In Sec.|IT] 
we will give a brief overview of the SIDIS production 
of one and two hadrons including the Sivers effect. In 


Sec. Ill we present our predictions for one hadron Sivers 
SSAs for both CLAS12 and EIC experiments, while in 
Sec. |IV| we present the similar results for the two hadron 
production. Finally, in Sec. [V] we present our conclusions 
and outlook. 


II. THE SIVERS EFFECT IN SIDIS 


In the SIDIS process a high energy lepton with momen¬ 
tum l, scatters off a nucleon target, N , producing several 
final state particles along with the scattered lepton with 
momentum V 


£(l) + N(P N ,S)^£(l') + X, (1) 

where P,y and S denote the momentum and the spin four- 
vectors of the nucleon. The leading-order DIS process in 
the factorized framework, where the lepton hard scatters 
of a single parton inside of a nucleon via single virtual 
photon exchange 7 *, is schematically depicted in Fig. [I] 
Here the momenta and the spin four-vectors of the initial 
and the scattered partons are denoted by k, s and k', s'. 

In describing the relevant kinematic variables and 
cross-section, we adopt the 7 * — N center of mass frame 
depicted in Fig [2] The z axis is taken along the direction 
of the virtual photon three-momentum < 7 , and the x axis 
is along the component of the lepton’s momentum trans¬ 
verse to q. The transverse components of the momenta 
in this frame are defined with respect to the z axis with 
subscript t and the transverse momenta with respect to 
the fragmenting quark’s direction with subscript j_. 



FIG. 1. The leading order diagram for one hadron production 
in the current fragmentation region of SIDIS. 
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FIG. 2. 7 * — TV center of mass frame. 


The usual DIS variables are defined as: 

q = i-r, Q 2 = -q\ X = ^S—’ ( 2 ) 

2 Pn ■ q 

y =^ L1 f W 2 = (P N + q) 2 , 

fN ' l 

where Pn is the momentum of the nucleon and x is the 
Bjorken variable. It is also straightforward to see that in 
the 7 * — N system the transverse momenta of the initial 
and fragmenting quarks are the same, k' T = kr, and 
qr = 0, by definition. 

The Sivers effect describes the modification of the 
TMD PDF of unpolarized partons inside a transversely 
polarized nucleon, fS(x, hr), with respect to that for the 
unpolarized nucleon ff(x,kr) (we omit the Q 2 depen¬ 
dence of these PDFs for notational convenience). This 
modification is quantified by the Sivers PDF, fff!(x, hr), 
which arises from the correlation between the active 
quark transverse momentum, kr, and the transverse po¬ 
larization of the nucleon, St- The polarized PDF is then 
expressed as (see e.g., mm) 

f*(x, kr) = fi(x, k T ) + kr )’ ( 3 ) 

where M denotes the nucleon mass and the subscript 3 
denotes the z component of the vector. 

A. The Sivers effect in single hadron production 


In the single hadron production in SIDIS, 

e(!)+N(P N ,S)-¥t(l') + h(P)+X, ( 4 ) 

one of the final state hadrons, h, with momentum, P , 
is detected along with the scattered lepton, Z, carrying 
momentum V, as schematically depicted in Fig. [l] Here 
two of the frequently used scaling variables are 


z = 


Pn-P 
Pn ■ q ’ 


(5) 


the fraction of the maximum possible longitudinal mo¬ 
mentum of the produced hadron in the 7 * — N center- 
of-mass CM system (Feynman x). In QCD factorized 
approach in the so-called current fragmentation region 
(CFR) defined as xf > 0, this hadron is primarily pro¬ 
duced by the fragmentation of the scattered parton. Al¬ 
though, within a dynamic string fragmentation picture, it 
is hard to make a strict distinction between the hadrons 
produced by the scattered parton’s fragmentation and 
those by the nucleon remnant’s fragmentation, especially 
for hadrons with relatively small momentum components 
along the lepton-nucleon axis in the 7 * — N CM system. 
Nevertheless, a formal distinction is assigned by consid¬ 
ering hadrons with xf > 0 in the CFR and those with 
xf < 0 in the target fragmentation region (TFR). 

Then, in the leading order approximation of QCD fac¬ 
torized approach for hadrons produced in the CFR, the 
relevant cross section can be expressed as [ 2 T] 


dx dQ 2 d^ s dz d 2 P T C{ ■ X, + <t S ) , 

4 = E e « / d 2 k T flD h q , 

q J 

°s = E e «/ d2kT ^W~ sin ^ fc ~ Vs)ftr D 1q, 


C(x, Q 2 ) 


a 2 (l + (l ~y) 2 ) 

Q A 


(7) 

( 8 ) 

(9) 

( 10 ) 


where ipu and ps are the azimuthal angles of the quark’s 
transverse momentum, kr, S the nucleon spin and a 
the fine-structure constant. D^ q (z, P ]_) denote the TMD 
fragmentation functions (FF), where P± is the transverse 
momentum of the hadron with respect to the fragmenting 
quark, acquired in the hadronization process. 

The Sivers SSA for this process can be expressed as 


A h = r, f di Psdph Off sm(<ph - p s ) 
Slv f dpsdp h cr(j 


( 11 ) 


Note, that the Sivers SSA of Eq. is calculated us¬ 
ing the cross-section formulae ( |8|9[ ) entail that we only 
use the leading order approximations (in twist expansion 
and a s expansion for the hard scattering, etc.). This ap¬ 
proach considers corrections from gluon radiation by the 
initial or scattered quark, or when the virtual photon 
interacts with the nucleon via vector-meson-dominance 
(VMD) type processes, to be small. These assumptions 
might not strictly hold for real world experiments, and 
we will explore this issue using MC generators for COM¬ 
PASS and EIC energies in the next section. 


B. The Sivers effect in two hadron production 


describing the fraction of the virtual photon’s energy car¬ 
ried by the produced hadron in the target rest frame, and 


Xf = 


2 P z 
W ’ 


( 6 ) 


Here we consider the SIDIS process, where along with 
the scattered lepton we also detect two of the produced 
hadrons hi and / 12 , with momenta P\ and P 2 

i(l) + N(P n , S ) -a 1(1') + Zii(Pi) + h 2 (P 2 ) + X. (12) 
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FIG. 3. The leading order diagram for two hadron production 
in the current fragmentation region of SIDIS. 


or ipT respectively: 
da hlh2 


d 2 P T RdR 


= C(x,Q 2 ) 


&u, o 


(19) 


Pt R 

+s T ( TF aT ’° + 1 2 m CTr ’ 1 ) sin ^ T _ vs) 


da 


h\ Ji2 


P T dP T d 2 R 


= C(x, Q 2 ) 


<?U,0 


( 20 ) 


+St 


) sin (^fl - <Ps) 


The leading order factorized schematic diagram for this 
process is depicted in Fig. [3j where we assume both 
hadrons are produced in CFR. 

In the experimental and phenomenological analysis of 
two hadron SIDIS, the transverse components of the two 
hadron momenta, Pit and Pit, are often replaced by 
their linear combinations 


P T = Pit + Pit, 

(13) 

R— 2 Tit — Pit) , 

(14) 


where ajj,i,aT,i and an,i are the zeroth (i = 0) and the 
first (i = 1) harmonics of the cos (n(tpT — <Pr)) Fourier 
expansions of the corresponding structure functions. The 
superscripts h\h ,2 are omitted. 

Then the relevant Sivers SSA can be obtained for 
the sine modulation of the cross section with respect to 
<Pt,r - <Ps 


A h!h 2 = J' d Vsdip a h s lh2 sin(y» - <p s ) 
StV f dipsdtp a'u lh2 


( 21 ) 


III. MONTE CARLO PREDICTIONS FOR ONE 
HADRON SSA 


where the corresponding azimuthal angles are denoted as 
<pt and ipn. The cross section terms relevant to the Sivers 
effect in this process can then be expressed as mm 


da hlfl2 


= C(x,Q 2 )(a^ h2 +a h s ^), 


dx dQ 2 dip s dzi dzi d 2 Pt d 2 R 


Jl}h2 =Y. 4 / d2k T fi 


Y h2 = 


2 f d 2 k T ^ T ! fcr] 3 /iT 9 Dhl 


M 


IT lq 


(15) 

(16) 

(17) 


where the virtual photon energy fractions z\ and zi car¬ 
ried by two hadrons are defined similarly to the one 
hadron case. D^ h2 is the fully unintegrated unpolar¬ 
ized dihadron fragmentation function. The terms ajj 
and as depend on x, Q , z\, Z 2 , Pt, R and Pt • R = 
PtRcos(ipt ~ <Pr)- The Sivers term can be expressed as 


as = St 


Pt . , . R . . . 

CTT 17 sm ^ T _ Vs) + CTfl — sm(^ fl - ips) 


(18) 


The cross section in Eq. (15) is usually measured after 
integrating over the azimuthal angle (and often magni¬ 
tude) of the relative or total transverse momentum, ipR 


The Monte Carlo (MC) event generators are im¬ 
mensely valuable tools for both phenomenology and ex¬ 
periment in hadronic and high energy physics. They al¬ 
low one to understand various features in the complicated 
DIS processes that are induced by the strong nonpertur- 
bative as well as perturbative effects, kinematical limita¬ 
tions coming from both underlying process and detector 
performance used in a particular experiment. Another 
advantage provided by MC event generators, versus sim¬ 
ply employing the analytical expressions for the relevant 
cross sections, is the full coverage of the available phase 
space for produced particles. Furthermore, the dynami¬ 
cal full hadronization models employed in MC generators 
provide a unified description of the process in all regions 
of xf, where the formalism for SIDIS in CFR for xf > 0 
involves different functions (PDFs and FFs) to those in 
TFR for xf < 0 (Fracture Functions [25M27] h 

Here we start with the well-established MC generator 
PYTHIA6 and modify it (mPYTHIA) to include the Sivers 
effect, which yields a modulation of the azimuthal angle 
of the struck quark’s transverse momentum with respect 
to the nucleon’s spin. For this we use the phenomeno¬ 
logical parametrization of the Sivers PDF of Refs. UM 
extracted from fits to SIDIS SSAs measured in the HER¬ 
MES 0 ES] and COMPASS [29} [301 experiments. The 
modified version of the LEPTO MC generator, mLEPTO, 
has already been used in previous work [HHninnsg 
to estimate the SSAs in both single and dihadron SIDIS 
production, mostly in the COMPASS kinematics. For 
the present study PYTHIA6 has been chosen as a more 
modern generator that among other new features uses 
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double precision to operate with the particle momenta 
which is better suited for EIC studies. 

The Sivers effect has been implemented in mPYTHIA 
by modifying the part of the PYTHIA6 code that gen¬ 
erates the azimuthal angle of the struck quark in the 
7 * — N CM system, where the x and kx have already 
been selected according to the cross-section formula for 
the lepton-nucleon DIS scattering and the selected in¬ 
trinsic transverse momentum parametrization (we chose 
Gaussian), respectively. Here, instead of the uniform dis¬ 
tribution for the case of an unpolarized nucleon, we sam¬ 
ple the azimuthal angle according to the Eq. ([3]). We use 
one of the standard parametrizations for the unpolarized 
PDF accessible in PYTHIA6. We take the Sivers function 
parametrization from Refs. [S, It)] (’’DGLAP evolution” 
set) and slightly modify it to best reproduce COMPASS 
measurements of the corresponding asymmetries. The 
latter is expressed as a product of the unpolarized PDF 
and two functions, where one depends only on x and the 
other only on kx- The parametrization actually involves 
only these two multiplicative functions. Further, we use 
a Lorentz transform to determine Sx in the 7 * — N CM 
system in order to calculate the Sivers term. These mod¬ 
ulations of the struck quark azimuthal angle are then 
transferred to the final state hadrons by the hadroniza- 
tion process. Moreover, these azimuthal modulations are 
also transferred to the nucleon remnant via momentum 
conservation. Note that here we only study the Sivers 
effect for the light quarks, as currently there are no avail¬ 
able parametrizations for the Sivers PDFs of a gluon or 
for heavy flavor quarks. 

Here we use mPYTHIA to study the Sivers SSAs for sin¬ 
gle hadron production in SIDIS at the soon to be oper¬ 
ational CLAS12 and planned EIC experiments. In this 
initial study, we do not include the TMD evolution of the 
Sivers PDF, as it would require major rewrites of these 
MC codes to also include the TMD evolution in the un¬ 
polarized PDFs and the hadronization model. Moreover, 
here we only use some basic kinematics limits coming 
from the angular coverage of the detectors in each ex¬ 
periment, without taking into account the realistic ac¬ 
ceptance effects in these experiments. Nonetheless, one 
should keep in mind that our results are aimed at estab¬ 
lishing order-of-magnitude estimates for the observables. 

First we performed a cross-check between mPYTHIA and 
mLEPTO for the COMPASS kinematics, where a trans¬ 
versely polarized proton asymmetry has been measured 
using Eft = 160 GeV muon beam [5]. In MC simula¬ 
tions we allowed only the leading order DIS process and 
excluded any parton showering effects. The results of 
this simulation, along with the COMPASS data [5] and 
mLEPTO results are depicted in Fig. [4] We see that there 
is overall good agreement between mLEPTO and mPYTHIA, 
with both describing the experimental data reasonably 
well. The small differences arising in ^-dependence be¬ 
tween the two MC generators are due to the differences 
in the hadronization modeling and parameters. 



X 


(a) 




/', (GeV) 

(c) 

FIG. 4. Sivers SSAs vs (a) x, (b) y and (c) Pt from COM¬ 
PASS data, and calculated with mLEPTO and mPYTHIA. 

A. The experimental kinematics 

We use the kinematic limits and cuts for CLAS12 de¬ 
scribed in Ref. [23], where an 11 GeV electron beam scat¬ 
ters from transversely polarized proton target. The rel¬ 
evant ranges for the kinematical variables are 1 GeV 2 < 
Q 2 < 6.3 GeV 2 and 0.075 < x < 0.532. The polar angle 
(in the lab frame) of the scattered electron 7° < 6 e ' < 35° 
and the polar angles of the hadrons 7° < 9^ < 100° are 
imposed by the detector coverage. We additionally ap¬ 
ply a cut W 2 > 4 GeV 2 to select mostly DIS events. For 
single hadron SSAs, we impose 0.2 < z < 0.8 . For di¬ 
hadron pairs we impose 0.4 < z < 0.8 for the total energy 
fraction of the pair. Additionally, to exclude the exclu¬ 
sive baryon resonance production region, we require the 
missing mass for both one- and two-hadron production, 
M( e +p)-(e'+h) and M {e+p) _ {e , +hl+h ^, to be larger than 
1.5 GeV. 
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The EIC [2D] proposal has several possible implemen¬ 
tations at different existing and planned facilities. Here 
we choose the kinematics of the experiment based on the 
RHIC facility [M], where electron and polarized proton 
beams collide. The predictions for EIC have been carried 
out for two sets of electron l and proton Py beam mo¬ 
menta, thus we will denote the corresponding energies in 
the text below as l x P.y and skipping the units always 
taken to be in GeV (for example 5 x 50). The following 
kinematic cuts are used: 

0.01 < y < 0.95, Q 2 > lGeV 2 , W > 5GeV, (22) 

as well as those described in Ref. [34] for SIDIS. Addi¬ 
tionally, we identify the hadrons in the CFR as those 
with xf > 0 in the 7 * — N CM system. 

B. EIC toy model studies 

We employ a toy model for the Sivers PDF to study 
the general features of the generated SSA. In particu¬ 
lar, it is interesting to study the relative contributions 
of the sea and valence quarks to Sivers SSAs in vari¬ 
ous regions of x and z. This will help to identify the 
most relevant experimental regimes for studying sea or 
valence Sivers PDF. For a toy model Sivers PDF, we 
choose it proportional to the unpolarized PDF of the 
corresponding quark flavor, with a large coefficient to 
maximize the generated SSA while satisfying the posi¬ 
tivity of the polarized PDF. The plots in Fig. [5] depict 
the results for 5 x 50 SSAs’ dependence on (a) x and (b) 
z. Here the label ’’Full” denotes the calculations where 
we set the Sivers terms using the toy model expression 
(x, kx) = 0.9 fi(x, kx) for all the light fla¬ 
vored quarks, while for calculations with the label ”Val” 
we only use the toy model Sivers terms for the valence 
u and d quarks and set them to zero for the sea quarks. 
Finally, the label ”Sea” denotes the scenario where we 
use the toy model Sivers terms for the sea quarks and 
set them to zero for the valence quarks. It is apparent 
that the valence quarks produce the predominant part of 
the signal in the large- 2 ; region, while the sea does in the 
small -x region, as one would naturally expect. Also, for 
all the cases, the SSAs increase as a function of z, as one 
naturally expects that the hadrons with a large fraction 
of the fragmenting quark’s energy also carry a large frac¬ 
tion of its transverse momentum. Thus their azimuthal 
modulations become larger. Finally, the plots labeled 
” Val+Sea”, indicating the simple sum of the results when 
we have sea-only and valence-only contributions. These 
coincide with the ’’Full” simulations, serving as a simple 
cross-check. 

Next, we study the dependence of the valence and sea 
quark contributions to SSAs for several different energies, 
as depicted in Fig. [g] The results where we set all the 
light quarks’ Sivers PDF to be the model value show 
little energy dependence, as only the increase in charm 
sea quark PDFs at higher CM energies slightly reduces 



(a) 



(b) 


FIG. 5. EIC toy model SSAs vs (a) x and (b) 2 for 5 x 
50 SIDIS kinematics for 7 r + . Comparison of sea and valence 
quark contributions in pure DIS events. 

the corresponding SSAs (here we do not allow a Sivers 
effect for the charm quarks). The valence and sea quarks, 
as expected, trade strength as the CM energy increases, 
and naturally the two results 5 x 250 and 30 x 50 with 
very similar CM energies also produce very close results 
for the corresponding SSAs. Thus one can conclude, that 
the high energy collision runs at EIC would be much 
more sensitive to sea quark (and by the same to the logic 
gluon) Sivers PDF, while the relatively low energy ones 
will most effectively probe the valence quark region. 

Studying the Sivers asymmetry in SIDIS in the TFR, 
xp < 0 , provides yet another avenue for accessing the 
spin-orbit correlations in the nucleon via so-called frac¬ 
ture functions H3I23- Note, that in the leading order 
QCD factorized approach, based on the fracture func¬ 
tions formalism, the only single spin asymmetry allowed 
in the TFR is the Sivers-like one. Thus, by measuring 
both this and the Collins-like [sin^^,+ <ps)] asymmetries 
(the latter predicted to be vanishing) one can probe fac¬ 
torization for SIDIS in TFR. Within the string model for 
hadronization adopted in the MC event generators, the 
inclusively measured hadron can be considered as orig¬ 
inating from the target remnant side of the string. If 
the struck quark’s transverse momentum is modulated 
by the Sivers effect, then so should be the transverse mo¬ 
mentum of the remnant with a relative phase 7 r via mo¬ 
mentum conservation. These modulations will be then 
transferred to the transverse momenta of the hadrons 
produced by the remnant, manifesting in the correspond- 
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(a) 



(b) 



(c) 


FIG. 6. EIC toy model SSAs for n + for a range of beam 
energies. Comparison of (a) the full, (b) valence-only and 
(c) sea-only Sivers effect induced SSAs are shown for various 
energies vs z, where we only allow for the DIS process. 


ing SSAs. Moreover, measuring the average transverse 
momentum square of the hadrons in the negative re¬ 
gion, xf < —0.2, allows greater sensitivity to accessing 
the quark average transverse momentum than in region 
xp > 0.2 [35]. The plot in Fig. [7] depicts the Sivers SSA 
for charged pions and kaons as a function of xf ■ We no¬ 
tice that the Sivers SSA is also nonvanishing in the target 
fragmentation region. Moreover, it changes sign when xf 
crosses zero, as expected from the arguments about the 
relative signs of the target remnant and struck quark’s 
transverse momenta. Thus we have demonstrated that 
the Sivers SSA measurements in TFR are as viable as 
in CFR, modulo any specific experimental setup limiting 
the hadron detection at negative Xf- 
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FIG. 7. EIC toy model SSAs for 5 x 50 SIDIS kinematics 
for charged pions and kaons vs xf- The Sivers asymmetry is 
present both in the current and target fragmentation regions. 

C. The effects of non-DIS processes and parton 
showering on A 8 ™ 

In this article we adopt the leading order approxi¬ 
mations for calculating the Sivers SSAs predictions for 
CLAS12 and EIC, where we ignore any possible non-DIS 
processes or gluon radiation corrections either by struck 
or fragmenting quarks (the initial- and final-state par- 
ton showering). Moreover, the same approach has been 
adopted when fitting the experimentally measured SSAs 
to extract the Sivers PDF [5; 'TDJ [35]. Thus it is very 
important to at very least make a qualitative study of 
the accuracy of this leading order treatment. A rigor¬ 
ous treatment of this problem would require that one 
include all the terms in the cross section for all the rel¬ 
evant processes, which would require the calculation of 
both unpolarized and polarized QCD and QED radia¬ 
tive corrections to DIS process to several orders in the 
corresponding couplings and take into account possible 
interference. Such an approach will additionally require 
extensive modifications of the existing MC generators to 
incorporate all of the new physics. Here we choose to use 
the less strict but more practical approach of using the 
PYTHIA6 generator with its built-in capabilities for treat¬ 
ment of the above effects to test the applicability of the 
leading order approximation. In PYTHIA6, the polarized 
parton showering is not included, while the unpolarized 
showering is assumed to occur either before or after the 
hard photon-parton scattering, with no interference ef¬ 
fects. Nevertheless, we think that a first step needs to 
be taken and the current machinery would suffice for an 
initial estimate. 

We start with the COMPASS kinematics and study 
four different cases. First we allow only the DIS single 
photon exchange process. Then, we also allow either par- 
ton showering (labeled ”+Show”) or non-DIS processes 
(labeled ”+NonDIS”) such as VMD and resolved pho¬ 
ton contributions. Finally, we allow both parton show¬ 
ering and non-DIS processes alongside the DIS process. 
The details of the particular methods used in PYTHIA6 to 
model all these processes can be found in Ref. m- 

The plots in Fig. [8] depict the x, z and Pp dependence 
of the Sivers asymmetry as we include various contribu- 
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FIG. 8. COMPASS SSAs for vs (a) x, (b) 2 and (c) 
Pt when various contributions beyond the DIS process are 
included. 


FIG. 9. EIC(5 x 50) SSAs for 7r + vs (a)x, (b) z and (c) 
Pt when various contributions beyond the DIS process are 
included. 


tions. Also depicted are the COMPASS measurements 
and the results obtained from the mLEPTO generator. Here 
we see that the non-DIS and parton showering contribu¬ 
tions significantly decrease both the x and z dependence 
of the SSAs, as the denominators of the SSAs, Eq. ©> 
with the ’’unpolarized” terms increase. It is also curi¬ 
ous, that the parton showering effects shift the strength 
of the SSAs from the low to high Pt region. Intuitively, 
this can be understood as the fragmenting quark being 
given an additional kick, uniform in </?, in the transverse 
direction by emitting a gluon or a photon, thus increasing 
the average Pt of the produced hadrons. Nevertheless, 
all of the results describe the current data reasonably 
well since it has significant uncertainties. On the other 
hand, there might be a potential impact from these find¬ 
ings on the phenomenological extractions of the Sivers 


PDF 0 EH 13B] from experimentally measured SSAs by 
HERMES and COMPASS. One can readily see from the 
Eq. (35) in [TO] , that the LO expression for the Sivers SSA 
is used in the fits would underestimate the relative con¬ 
tribution of the Sivers term in the cross section. Thus we 
deduce that the current extractions of the Sivers PDFs 
are most likely to significantly underestimate the true 
values because of the effects described. The reevaluation 
of these parametrizations might pose a significant chal¬ 
lenge, similar to those described earlier in this subsection. 
On the other hand, here our goal is to predict the SSAs 
to be measured at the future experiments. Thus we need 
to first understand what is the relative contribution of 
these processes at the kinematics of interest. 

The plots in Fig. [9] are the analogues of those in Fig. [8] 
calculated for EIC 5 x 50 kinematics. Here we see trends 
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FIG. 10. Ratios of full to DIS-only results for rates ( N ) and 
SSAs (A Si ”) for EIC(5 x 50) and COMPASS vs (a) x, (b) a 
and (c) Pt • 


in the behavior of the SSAs similar to those in COMPASS 
kinematics. This hints that the relative changes of the 
SSAs might be very similar at both COMPASS and EIC 
kinematics. To test this hypothesis we have calculated 
the ratios of both hadron count rates and SSAs calculated 
including all the non-DIS and parton showering processes 
to those calculated only including the DIS process. 

The plots in Fig. [To] depict the ratios of full and DIS- 
only results vs (a) x, (b) z and (c) Pt for EIC ( 7r + at 
5 x 50) and COMPASS ( h + ). We can see striking simi¬ 
larities for both COMPASS and EIC results, if we do not 
take into account the small-a; region for the SSAs, where 
in the ratios of two very small numbers the small statis¬ 
tical uncertainties create a large spread of results. These 
results are encouraging, as we can still use the Sivers 
PDF parametrization m to obtain good estimates of 


the SSAs using only the DIS process. The reason is that, 
to obtain the COMPASS SSA induced by the DIS-only 
contributions to both unpolarised and Sivers terms, the 
measured SSA should be divided by the ratio of the full to 
DIS contributions for COMPASS depicted in Fig. 10 On 
the other hand, to predict the size of the SSAs at EIC, 
the estimates that include only the Sivers PDF contri¬ 
butions should be multiplied by the ratio of the full to 
DIS contributions for EIC, also depicted in Fig. 10 We 
can readily see from Fig. [lO] that these two factors cancel 
each other to a very approximation. Moreover, our ini¬ 
tial calculations [TS] that only included the leading order 
DIS process were reproducing the measured COMPASS 
data well. A rigorous study of the effect of polarized and 
unpolarized parton showering, as well as TMD evolution 
of the Sivers PDF in MC simulations for calculating the 
observable SSAs is left for future work. 


D. EIC projections 

Here we present the results for the EIC simulations 
using the kinematical cuts given in the beginning of the 
section. In simulations we only allow pure DIS scattering. 
We present results for two sets of energies of the electron 
and proton beams, 5 x 50 and 5 x 250, that provide a 
reasonable variation in the CM energies. 

The plots in Fig. [TT] depict the SSAs for charged pseu¬ 
doscalars as functions of (a) x, (b) 2 and (c) Pt for the 
5 x 50 configuration. We see that the SSAs for K + are 
large, while those for tt~ and K~ are relatively small for 
all cases. This can be understood in terms of a signifi¬ 
cant u —> K + FF that has a bump at relatively large z 
(z ss 0.7), while u —> ir + drops sharply with increased 
x (though it stays larger in magnitude than that for the 
K + ). Thus, on average a larger portion of the fragment¬ 
ing quark’s transverse momentum is transferred to the 
K + , leading to the large SSA. Additionally, the contri¬ 
bution of the s quark to the K + asymmetry is quite 
large due to the parametrization of the Sivers function 
for s and s. which was driven by HERMES and COM¬ 
PASS data. The contribution of the strange quarks to tt + 
is negligible, as the corresponding fragmentation func¬ 
tions are very small compared to those for the up quark. 
We have checked explicitly, that when the strange Sivers 
PDF is set to zero, only the K + SSA is affected, it is 
reduced and becomes almost equal to that for the ir + . 

Next we want to test how these predictions would vary 
for increased CM energy, corresponding to the 5 x 250 
configuration of beam energies. The plots in Fig. [12] de¬ 
pict the results for the SSAs for ir + as functions of (a) x, 
(b) 2 and (c) Pt for this energy. We see that the SSAs 
in general decrease as we increase the CM energy. This 
is because at larger energies we sample lower x values, 
where the sea quarks give the largest contributions to 
the SSAs. Note that we have used the parametrizations 
extracted from COMPASS and HERMES data, where 
the Sivers PDF for the sea is not very well constrained 
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(a) 



(b) 



FIG. 11. EIC SSAs vs (a) x, (b) z and (c) P T for 5 x 50 SIDIS 
production of charged pseudoscalars. 

and comes out small when compared to valence. Thus we 
can only conclude from these plots that if the sea Sivers 
PDFs are small compared to the valence ones, then the 
Sivers SSAs should become smaller and smaller as we in¬ 
crease the energy. On the other hand, measurements at 
higher energies would help to probe the small x region 
and better measure the sea Sivers PDFs. 

The plots in Fig. [13] depict the results of the SSAs 
for charged pions and kaons versus x f, analogous to 
those depicted in Fig. [7] for the toy model. We see 
that the nonzero Sivers SSAs are also present in the 
TFR (xp < 0), with opposite sign to those in the CFR 
(xf > 0) except for n~. This and the nontrivial shapes 
and unequal sizes of the SSAs for the other mesons are 
governed by the Sivers PDF parametrizations and the 
complicated hadronization dynamics in the CFR and 
TFR regions. Thus, measuring these SSAs over the en¬ 



(a) 



(b) 



FIG. 12. EIC SSAs vs (a) x, (b) z and (c) P T for 5 x 250 
SIDIS production of charged pseudoscalars. 

tire range of Xf (or rapidity) will provide additional in¬ 
formation to improve our understanding of the complete 
hadronization of both the knocked out parton and the 
target remnant. 

One of the important challenges in phenomenological 
fits is to disentangle the flavor dependence of the Sivers 
PDF from the measurements for a small set of final par¬ 
ticle types. The difficulty is that in the corresponding 
cross section the Sivers PDF is convoluted with the frag¬ 
mentation function to the specific hadron and these are 
summed over the flavor of the parton. Thus having SSA 
measurements for as large a number of final hadron types 
as possible would allow us to remove some or all of the 
assumptions about the flavor dependence of the corre¬ 
sponding Sivers PDFs. In experiment though, the prac¬ 
tical limitations of the detectors limit this ability to only 
a few hadrons, such as charged pions and kaons. Thus 
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FIG. 13. EIC model SSAs for 5 x 50 SIDIS kinematics for 
charged pions and kaons vs xf ■ The Sivers asymmetry is 
present both in the current and target fragmentation regions. 

it is important to gain phenomenological understanding 
of the relative contributions of each of the parton flavors 
to the Sivers term in the cross section. In Fig. [14] we 
depict the results for the rates of struck quark flavor for 
the selected events that contribute to the count rates of 
a given hadron for charge pion and kaon final states and 
the two different sets of energies. We see that the charm 
sea becomes significant at large CM energies, though we 
set the charm Sivers function to be zero. Thus this should 
be in part responsible for weakening of the SSA signals 
with increasing energies, as we saw earlier. We have seen 



(a) 



(b) 

FIG. 14. EIC rates for the struck quark flavor for the selected 
events in SIDIS productions of charged pions and kaons for 
(a) 5 x 50 and (b) 5 x 250 energies. 



(a) 



(b) 



FIG. 15. Predictions for the Sivers asymmetry in charged 
hadron production at CLAS12 versus (a) x, (b) z and (c) Pr¬ 
ill the toy model calculations that with the increased en¬ 
ergy there was a slight reduction of the SSAs, even when 
all the light quarks were assigned a large Sivers PDF. We 
should keep in mind here that these integrated rates hide 
a much more complicated picture when considering the 
full x, z and Pt dependences. 


E. Predictions for CLAS12 

In this subsection we study the single hadron SSAs 
for the upcoming CLAS12 experiment Ezl> again using 
the mPYTHIA MC event generator with the parameters 
slightly tuned to best suit the kinematics for the CLAS12 
experiment at 11 GeV, that has been described in the 
beginning of this section. 

The mPYTHIA result for the one hadron Sivers SSAs for 
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perimental uncertainties. Further, development and test¬ 
ing on real world data of a complete analysis tool chain 
for extracting the Sivers PDF, that accounts for all the 
non-DIS and parton showering effects for CLAS12, would 
be of high importance both for CLAS12 and future EIC 
experiments. 


IV. MONTE CARLO PREDICTIONS FOR TWO 
HADRON SSA 


FIG. 16. Predictions for the relative rates for the flavor of the 
struck quark that produces charged hadrons in MC events at 
CLAS12, with all the relevant kinematical cuts. 
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FIG. 17. Predictions for SSAs for charged pions and kaons vs 
xf at CLAS12. The Sivers asymmetry is present both in the 
current and target fragmentation regions. 

charged pions and kaons as functions of x, z and Pt are 
presented in Fig. 15 (a), (b) and (c), respectively. In 


Fig. 16 (drawn in the log scale) we show the analogous 


plots to Fig. [14] with the rates of the produced hadrons 
after all the relevant cuts have been applied. We note 
that the Sivers SSA for 7r + is large and positive, but 
the SSA for K+ is even bigger. On the other hand, the 
rates for K + production are about an order of magnitude 
smaller than for 7r + . Thus much higher statistics would 
be required to measure K + SSAs versus those for . 
The SSAs for tt~ and K~ are relatively smaller, with 
the latter being produced at a rate that is much less 
than even K + , making it very challenging to measure at 
CLAS12. 


Finally, the plots in Fig. [17] depict the results of the 
SSAs for charged pions and kaons versus xf ■ It is inter¬ 
esting to note that the SSAs here are larger in magnitude 
than those plotted for EIC in Fig.[T3j thus they should be 
easier to access in the experiment. Additionally, we no¬ 
tice here that the kinematical cuts employed for CLAS12 
(on the missing mass, etc) limit the xp region of these 
mesons to about ±0.75. 

Here we would underline the importance of the new, 
higher-statics measurements of the Sivers effect in one 
hadron SIDIS. This will allow us to improve the extrac¬ 
tions of Sivers PDFs, especially in the valence region, as 
well as allowing us to shrink the error bars of the corre¬ 
sponding parametrizations coming from the current ex- 


The one hadron, SIDIS production off a transversely 
polarized target has been long established as one of 
the principal processes for measuring the Sivers PDF. 
Nonetheless, the phenomenological extractions of the 
Sivers PDF still remain challenging. This is especially 
true for the flavor dependence, where we can only mea¬ 
sure a limited number of different hadron types in the fi¬ 
nal state. Two hadron SIDIS production [14UI6] provides 
a complimentary approach to extract the Sivers function 
with very different systematic errors both from theory 
and experiment. Here again the relevant term of the 
cross section is a convolution of the Sivers PDF with an 
unpolarized dihadron fragmentation function, summed 
over the type of the parton. This approach is promising, 
because with two detected hadrons in the final state we 
have a much wider basis for extracting the flavor depen¬ 
dence of the Sivers PDF. For a given number of types 
of the detected hadrons Ay,, the number of possible di¬ 
hadron pairs is given by • (Nh + l)/2, as we can have 
pairs with the same types of hadrons. Thus it would be 
important to estimate the viability of such dihadron SSA 
measurements at both CLAS12 and EIC using our MC 
generators. Here we can measure SSAs corresponding 
to modulations with respect to VSiv,T = Vt ~ Vs and 
VSiv,R = VR ~ Vs [see Eqs. ( 19|20 1], that we will simply 
refer to as ”T” and ”R” modest 


A. Predictions for CLAS12 


1. Dihadron SSAs 


Here we consider the dihadron pairs of pions and kaons. 
In identical pairs (e.g. 7r + 7r + ), we choose hi as the hadron 
with the larger value of Zi in the pair. Figures. [j~8| show 
the results for both ”T” and ”R” type SSAs for dihadron 
pairs as functions of x. We see that the pairs involving 
K + have the largest SSAs. Further, the modulations in¬ 
volving the vsiv,T a re larger in magnitude than those for 
the same pair involving VSiv.R- Both SSAs, but espe¬ 
cially those for vsiv,R , can be enhanced significantly be 
choosing asymmetric cuts on the momenta of the pair, as 
has been already shown in COMPASS kinematics in [T5] 
and for EIC kinematics later in this sections. The de¬ 
pendences of these SSAs on the total energy fraction 2 
and invariant mass M of the pairs are presented in 
Figs. 19 and 20 It is interesting to note that for the 
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FIG. 18. Predictions for Sivers SSA for charged (a) pion 
and (b) kaon-inclusive pair production off proton target at 
CLAS12 versus x. 

same charged pion pairs, the ”T” and ”R” SSAs trade 
strength as we increase h 2 ■ 

Again, we can learn about the relative rates of the dif¬ 
ferent hadron pair production and their dependence on 
the struck quark’s flavor from the plots in Fig. [2TJ that 
are the analogues of those for the single hadron produc¬ 


tion in Fig. 16 tlere we note that tt^tt pairs are pro¬ 
duced at the highest rate, while the rate for K + K~ pairs 
is about 100 times smaller. Thus measuring kaon pairs 
would be more challenging than the pion pairs, but would 
be very important in understanding the flavor structure 
of the Sivers PDF. 


B. Predictions for EIC 


Here we present the dihadron SSAs for EIC kinemat¬ 
ics. First, we study the x dependence of ”R” and ”T” 
type SSAs for various charged hadron paris for 5 x 50 
kinematics, as depicted in Fig. 22 (a) and (b) respec¬ 
tively. Here again the SSAs involving the K + are the 
largest ones. Moreover, we see the different behaviors of 
these SSAs that would be very important in extracting 
the Sivers PDF from them. 


In Fig. 23 (a) and (b) we present the 5 x 250 results 
for both ”R” and ”T” type SSAs as functions of x for 
the same pairs as in Fig. [22] Again, similar to the one 
hadron case, we see decrease of the SSAs as we increase 



(a) 



z 


(b) 

FIG. 19. Predictions for Sivers SSA for charged (a) pion 
and (b) kaon-inclusive pair production off proton target at 
CLAS12 versus z. 

the CM energy. The asymmetric cuts on the momenta of 
hadron in the pairs should enhance these SSAs (especially 
in ”R” mode), as demonstrated for 7r + 7r^ pairs in the 
next subsection. 

The plots in Figs. [24] [25] depict the analogous results 
for both energies, but now as functions of z. Here the 
rise of most SSAs with z can be interpreted by simplis¬ 
tic arguments that the hadrons (and pairs) with larger 
z carry larger fractions of the struck quark’s transverse 
momentum. 

Finally, it would be important to once again examine 
the relative rates of the different hadron pair production 
and their dependence on the struck quark’s flavor, as de¬ 
picted on the plots in Fig. [26] Here we notice a significant 
increase in the rates of production from all the quarks 
compared to CLAS12 results shown in Fig. [2l] Thus 
EIC should provide a higher precision measurements for 
the relevant SSAs than CLAS12 for the same number of 
DIS events. This is of course due to the higher average 
multiplicities and energies of the produced hadrons. 


1. Dihadron SSAs with asymmetric cuts 

Following Refs. mug, here we also impose asymmet¬ 
ric cuts (zi > 0.3 and Pit > 0.3 GeV) on the momenta 
of the pair. These cuts enhance the SSAs, especially the 
”R” type, but reduce the rate of detecting such pairs 
and increasing the statistical noise. Thus the cuts might 
be optimized to gain the maximum signal to noise ratio 
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(a) 



(b) 


FIG. 20. Predictions for Sivers SSA for charged (a) pion 
and (b) kaon-inclusive pair production off proton target at 
CLAS12 versus Mh 1 h 2 - 



FIG. 22. EIC 5 x 50 results vs x for (a) ”R” and (b) ”T” type 
SSAs. 


for the particular experiment. Here we simply demon¬ 
strate the effects of the sample choice of the cuts on the 
dihadron SSAs. 


Here we present plots analogous to those from the pre¬ 
vious subsection, but now with the cuts applied. The 
dependences of the SSAs for various charged pairs on x 
and 2 are depicted in Fig. 27 (a) and (b), respectively. 
We note that in general, the SSAs are significantly en¬ 
hanced, especially those for ”R” type modulations. The 
tradeoff here is that the rates for producing the pairs are 
reduced by roughly three to four times. Thus a careful 
optimization for SSA enhancement/rate decrease can be 
made for particular experiment to achieve the best signal. 


( a ) 



(b) 

FIG. 21. Predictions for the relative rates for the flavor of the 
struck quark that produces charged (a) pion and (b) kaon- 
inclusive pair in MC events at CLAS12, with all the relevant 
kinematical cuts. 


V. CONCLUSIONS 

The measurement of the Sivers PDF remains one of 
the top priorities of the current and future experiments 
in hadronic physics, with SIDIS on a transversely polar¬ 
ized nucleon target being one of the most viable chan¬ 
nels. There are many challenges in achieving this goal, 
both in measuring the SSAs and extracting the PDF 
from them using phenomenological analysis of the data. 
Our goal was to use the state-of-the-art knowledge of 
the Sivers PDF to make predictions for the upcom¬ 
ing experiments at CLAS12 and the planned EIC. We 
used the well-established MC event generators LEPTO and 
PYTHIA6, modifying them to include the Sivers modula¬ 
tion of the struck quark’s transverse momentum using the 
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FIG. 23. EIC 5 x 250 results vs x for (a) ”R” and (b) ”T” 
type SSAs. 

phenomenologically extracted Sivers PDFs from HER¬ 
MES and COMPASS measurements. 

First we studied the single hadron SIDIS production in 
Sec. |III| We used a toy model to study the relevant con¬ 
tributions of sea and valence quarks to the SSAs for EIC 
kinematics at various energies. At the lower CM energies, 
the contributions of both sea and valence quarks were 
comparable in different regions of x, while the sea domi¬ 
nates the SSAs for the larger CM energies as one would 
intuitively expect. Next, we explored the contributions 
of non-DIS and parton showering effects on the SSAs, as 
in most of the phenomenological analyses these are omit¬ 
ted. We found a significant dilution of the Sivers SSAs 
when these processes were included in the MC simula¬ 
tions, strongly suggesting that the current extraction of 
the Sivers PDFs are significantly underestimated. Nev¬ 
ertheless, a comparative analysis of these effects both for 
COMPASS and EIC showed that these PDFs can be used 
to estimate the SSAs if we only include the DIS events, 
since the modifications of both unpolarized and polar¬ 
ized parts of the cross sections are very similar for both 
kinematics in this model. Then we presented the SSA 
results for various hadron production at EIC, and also 
explored their dependence on the CM energy. In the 
current parametrizations, the sea quark Sivers PDFs are 
small and poorly constrained, leading to decreasing SSAs 
with increasing CM energy. Similar results were also pre¬ 
sented for CLAS12. 

In Sec.[IV]we presented the SSA predictions for various 
charged hadron pairs in both EIC and CLAS12 kinemat¬ 


FIG. 24. EIC 5 x 50 results vs z for (a) ”R” and (b) ”T” type 
SSAs. 

ics. We have noticed that they are of a similar size to 
the single hadron ones, thus making them accessible ex¬ 
perimentally. Further, asymmetric cuts on the momenta 
of the hadrons in the pair allow one to significantly en¬ 
hance these asymmetries at the cost of reduced statistics. 
Moreover, the dihadron channel provides a large basis for 
disentangling the flavor dependence of the Sivers PDF. 
The measurement itself is complementary to the highly 
publicized two hadron studies of the transversity PDF via 
the interference fragmentation mechanism, with several 
current analyses of the experimental data and a num¬ 
ber of proposals for future experiments using CLAS12 
and SoLID at JLab. In fact, the dihadron Sivers analy¬ 
sis has been already included in the SoLID proposal [IS] 
for transversity measurements on a neutron target, that 
would provide very exciting test for the Sivers PDFs ex¬ 
tracted from proton and deuteron data. On the plreno- 
logical side, these dihadron measurements would comple¬ 
ment the single hadron measurements in the global data 
fits. The unknown, fully unintegrated dihadron fragmen¬ 
tation functions needed for such analysis also need study, 
but a very similar challenge exists also for the dihadron 
transversity measurements (there we need two dihadron 
FFs). Then these dihadron FFs can be studied in a sim¬ 
ilar fashion, such as using e + e _ data on dihadron multi¬ 
plicities, MC and model studies. 

The measurements of the SSAs for hadrons produced 
in the TFR will allow one to study the structure of 
the nucleon through the fracture functions. These ob¬ 
jects, though more complicated than the ordinary PDFs 
and FFs, will provide important new information. The 

















16 



(a) 



(b) 


FIG. 25. EIC 5 x 250 results vs 2 for (a) ”R” and (b) ”T” 
type SSAs. 

hadronization in mPYTHIA and mLEPTO is governed by the 
Lund string model, that differs from the usual QCD fac¬ 
torized approach that describes the hadron production 
in the CFR with a convolution of PDFs and FFs and 
in the TFR using two additional independent fracture 
functions. In particular, in mPYTHIA we take into account 
only the correlation between the nucleon’s transverse po¬ 
larization and the transverse momentum of the struck 
quark, which propagates to both CFR and TFR hadrons 
by the Lund fragmentation model. On the other hand, 
in the TFR factorized formalism there are two indepen¬ 
dent fracture functions describing the correlation of both 
the struck quark’s and the produced hadron’s transverse 
momenta with the transverse spin of the nucleon. Thus 
it would be crucial to test experimentally, whether the 
unified description of Sivers SSAs in both CFR and TFR 
can be achieved with mPYTHIA. Our results of the Sivers 
SSAs for charged hadrons showed a sizable signal in the 
TFR, comparable in size to that in the CFR. 

In conclusion, we have developed a tool based on very 
successful MC generators for studying the Sivers effect 
in experiments. Our Sivers SSA predictions for both 
CLAS12 and EIC indicate that these should be accessible 
experimentally and will provide information on the Sivers 
PDF in a wide range of kinematical variables. We noted 
in particular the power of dihadron studies in pinning 
down the dependence of the Sivers PDF on the flavor of 
the struck quark. We did not include the TMD evolution 
of the Sivers PDF that would weaken the signal at EIC 
energies. Further, a detailed analysis of the real world 




(b) 

FIG. 26. EIC 5 x 50 predictions for the relative rates for 
the flavor of the struck quark that produces charged (a) pion 
and (b) kaon-inclusive pair in MC events with all the relevant 
kinematical cuts. 

experimental limitations, such as detector acceptance ef¬ 
fects, statistics with the specific proposed time frame of 
data collection, etc, will of course create more challenges 
in their measurements but these considerations are out 
of scope of this article. We have also demonstrated using 
the mPYTHIA MC generator, that the current extractions 
of the Sivers PDF might be significantly underestimated. 
Thus improvements in estimating the true relative size of 
the Sivers term in the cross section would be needed in 
such future analyses. 
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(a) 



(b) 


FIG. 27. EIC results for Sivers asymmetry versus (a) x and 
(b) z with asymmetric cuts applied to the momenta of the 
hadrons in the pair. 
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